UV fixed-point structure of the three-dimensional Thirring model 
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We investigate the UV fixed-point structure of the three-dimensional Thirring model by means of the func- 
tional renormalization group (RG). We classify all possible 4-fermi interactions compatible with the present 
chiral and discrete symmetries and calculate the purely fermionic RG flow using a full basis of fermionic four- 
point functions in the point-like limit. Our results show that the UV complete (asymptotically safe) Thirring 
model lies in a two-dimensional coupling plane which reduces to the conventional Thirring coupling only in 
the strict large-A'f limit. In addition to the Thirring universality class which is characterized by one relevant 
direction (also at finite A'f), two further interacting fixed points occur which may define new universality classes 
of second-order phase transitions also involving parity-broken phases. The A'f-dependence of the Thirring fixed 
point sheds further light on the existence of an A'f-controlled quantum phase transition above which chiral 
symmetry remains unbroken for arbitrary large coupling, even though a definite answer will require a direct 
computation of competing orders. 



I. INTRODUCTION 

Three-dimensional relativistic fermionic systems have ex- 
tensively been investigated in the literature in a variety of 
scenarios. On the one hand, they are per se interesting field 
theories with unconventional features, on the other hand they 
allow for fascinating applications to condensed-matter sys- 
tems. In particular, three-dimensional quantum electrodynam- 
ics (QED3) or the Thirring model iH] are actively discussed, 
e.g., as effective theories describing different regions of the 
phase diagram of high-Tc cuprate superconductors [2-3] and, 
recently, the electronic properties of graphene jsHa]. Espe- 
cially graphene since its first synthesis in 2004 lUO'l is be- 
ing lively discussed in the rapidly growing literature on this 
subject, also because it offers the appealing opportunity of 
a comparatively simple experimental realization of some up 
to now unobserved quantum relativistic phenomena, such as 
the Klein paradox |11] or Zitterbewegung [12]; for reviews, 
see llT3i[l4i1 . More exotically, some features of 3(i relativistic 
fermion systems can serve as toy models, e.g., for the stan- 
dard model of particle phvsics il5[l. or a possible candidate 
for a "Theory of Everyt hing " [Tg]- 

However, OED3 lfl7l - l23ll and the three-dimensional Thir- 
ring model ll24l - [3lll are likewise intrinsically interesting, in 
particular because the ground states in these theories are ex- 
pected to show a sensitivity to the number of fermion flavors 
A^f. Several approximate solutions of the Dyson-Schwinger 
equations (DSE) predict that chiral symmetry breaking (^SB) 
is prohibited once A^f is larger than a critical value A^^' 
ll26l I281 i29i1 . A similar quantum critical phenomenon has also 
been identified in many-flavor QCD4 ll32l[33ll . being currently 
under intense investigation also because of its potential rele- 
vance for technicolor scenarios [34] and its implications for 
the thermal phase boundary ifssil . 

The search for the quantum critical point in the Thirring 
model has so far been rather challenging: different DSE ap- 
proaches, e.g., have yielded values between A^" ^ 3.24 12611 
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and Nf - 00 ||27[1 . Recent extensive lattice simulations now 
point to A^" ^ 6.6 lIsTil . For physical graphene and cuprates, 
the number of flavors is A^f = 2, such that the true value of 
A^" may be of profound importance for physical effects cor- 
responding to chiral symmetry breaking in the effective the- 
ories. In fact, a dynamically induced mass gap in the band 
structure of graphene, corresponding to a semimetal-Mott in- 
sulator quantum phase transition, could provide extraordinary 
applications for graphene-based electronics, offering a possi- 
ble new candidate material to take over from Si-based tech- 
nology 1 13]. 

In this work, we take a more fundamental viewpoint in or- 
der to investigate the amount of universality that can be at- 
tached to a possible value for A^^''. Naively, a universal answer 
for A^" may not be apparent as id fermion models are pertur- 
batively nonrenormalizable. Nevertheless, renormalizability 
of the 3(i Thirring model to any order in a large-A^f expansion 
has been shown in i24l-l27^. with a diagrammar being very 
similar to QED3. This has been taken as an indication that 
the Thirring model can indeed by defined nonperturbatively 
in 3d, providing the same amount of universality as any per- 
turbatively renormalizable theory. 

Universality and, more profoundly, UV completeness can, 
in fact, be analyzed within Weinberg's scenario of asymptotic 
safety [36., 371, in which a UV complete infinite-cutoff limit 
can be taken at a potentially non-Gau6ian fixed point. The re- 
sulting theories are predictive and exhibit universality, if the 
number of RG relevant directions at the fixed point are fi- 
nite. In this work, we re-examine renormalizability of the 3d 
Thirring model from this viewpoint by means of the functional 
renormalization group formulated in terms of the Wetterich 
equation l38ll . Using a purely fermionic RG as a first step, 
the UV fixed-point structure can indeed be mapped out, on 
the one hand confirming the large-A^f results, but also reveal- 
ing interesting deviations from the large-A^f asymptotics. For 
instance, we can identify a non-GauBian UV fixed point that 
defines a "Thirring universaUty class" which, however, corre- 
sponds to a pure Thirring coupling only in the strict large-A^f 
limit. In addition, we find further fixed points which may be 
associated with phase transitions (and corresponding new uni- 
versality classes), e.g., towards a parity broken phase. 
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This paper is organized as follows. In Sec. |ll] we analyze 
the discrete and continuous symmetries of the classical the- 
ory. We classify the different possible mass terms as well as 
all possible fermionic interaction channels in the point-like 
four-fermion limit. With the aid of the functional RG, we in- 
vestigate the purely fermionic RG flow in Sec. HUl followed 
by an analysis of the fixed-point structure of the given class of 
models for varying number of fermion flavors in Sec. |IVl In 
Sec. [Vl we classify the resulting RG trajectories to the long- 
range physics, pointing to the existence of various phases with 
dynamical mass generation, and obtain a first glance at a pos- 
sible mechanism for the formation of a critical flavor number, 
above which chiral symmetry breaking disappears. We give 
our conclusions in Sec.lVIl 

II. SYMMETRIES OF THE THIRRING MODEL 

The Lagrangian for the massless Thirring model in three 
EucUdean space-time dimensions is 

i: = ^«i#° + ^(^«r,'A''f , (1) 

satisfying Osterwalder-Schrader positivityQ The index a runs 
over A^f distinct fermion flavors. In three dimensions, we 
could use 2x2 matrices (e.g., the Pauli matrices) as an ir- 
reducible representation of the Euclidean Dirac algebra 

[rp,7v}^26fty, (2) 

but this representation does not permit a chiral symmetry. We 
therefore work exclusively with a 4 x 4 reducible representa- 
tion of the Dirac algebra; an explicit representation is given 
by 

r. = (i°^T| . = 1,2,3, (3) 

with {(Ti , 0-2, 1x3) being the 2x2 Pauli matrices. Thus tp" rep- 
resents a four-component Dirac spinor. Now there are two 
further 4x4 matrices which anticommute with all y^j as well 
as with each other, 

74 = 1^ 0) 75 = 71727374 = ^0 (4) 

The massless Lagrangian ([TJ then is invariant under the axial 
transformations 

U^,(l): i// ^ e'"^* iff, tfr^iffe"^', (5) 

Vy,(l): (A e'^''= i/r, ifr ^ ifie'^^\ (6) 



For Osterwalder-Schrader positivity | 39], we require invariance of the ac- 
tion under (generalized) complex conjugation defined by i//^ := iij/yj to- 
gether with reflection of the (euclidean) time coordinate (which we choose 
to be X3). For a detailed discussion of our chiral conventions, see e.g., 
Ref. HSl. 



as well as the vector transformations 

Ui ( 1 ) : i/r e"V, {f/^{f/ (7) 

( 1 ) : tfr^ e'^''-'= iff, ijj ^ ^ e"'^''« , (8) 

with :- vyAy^. For each flavor a - 1, . . . , A^f, the theory 
thus has a global U(2) symmetry with the Hermitian genera- 
tors Tj, j - 1 , . . . , 4, 

Tj - 1,74,75,745. (9) 

This symmetry transformations together with flavor rotations 
form a larger symmetry of the classical massless Lagrangian, 
corresponding to the group U(2A^f) with the (2A^f)^ generators 

Ai^Tj, i^\,...,N{^, = ..,4. (10) 

Here, {A\,. . . , Aj^i_{\ are the generalized A^f x A^f Gell-Mann 
matrices, and Aj^i :- In, is the identity. In other words, com- 
bining the A^f four-component spinors if/' - (ifr^, if/^)^ (each 
consisting of 2 two-component Weyl spinors i/'l, i^r) into one 
collective 4A^t -component spinor 

(consisting of 2A^f two-component Weyl spinors), the theory 
is invariant under 

^ W, ^ , ^( e U(2N{), (11) 

where the entries of the unitary (2A^f) x (2A^f) matrix 1/ are 
complex numbers times the 2x2 identity matrix I2. 

Because of the reducible representation of the Dirac alge- 
bra, discrete transformations can be implemented in various 
ways 126]. Let us define 

C : iff 'r^ (i/'Q)^, 1/' i-> -(C^i/')^ (charge conjugation), 
V : ifi <r-^ P^ip, ijj (^f I (parity inversion), 

7~ : i/' i-> T,jif/, if/ 1-> ifrT^ (time reversal), 

where C and !P are unitary and T is antiunitary. The argu- 
ments of the transformed fields are x :- (-xi, X2, xj.) in the 
case of parity, and x :- {xi,X2, -xj,) in the case of time rever- 
sal. Q, P(, and T,, are unitary 4x4 matrices given by 

Q = 5 [(1 + ^7274 + i(l - ^7275] , (12) 

= 5 [(1+^)7174 + i(l-Ori75], (13) 

7^,, = 5 [(1+^)71 +1(1-^)7273], (14) 

depending on the pure phases ^, rj with |^| = \(\ - \t]\ - 1. 
Recall that 71 and 73 are antisymmetric and purely imaginary, 
whereas 72, 74, and 75 are symmetric and real. We thus see 
that on the classical level our theory is invariant under any of 
the discrete transformations C, V, and T individually, irre- 
spective of the values of the phases ^, ^, and 77. 
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Let us now consider building blocks of the effective action, 
starting at the two-fermion level. There are, in fact, four pos- 
sible mass terms ij/Tj\l/ with Tj given in Eq. ^ and diagonal 
flavor structure. (We shall suppress the flavor index a as long 
as it is not needed.) However, the term i/'(im + m'y^)^ trans- 
forms under ^ i-> e'"'''' ij/ into a parity even mass term oc ^ij/ 
if a is chosen to satisfy 2a - arctan(m'/?M). The analogous 
statement holds for the mass term involving 75. More gener- 
ally, any mass term can be transformed by a U(2A^f) rotation 
intcH 



itj/ (m + my 45) iff. 



(15) 



A dynamically generated mass m Q spontaneously breaks 
the U(2A^f) symmetry down to a residual Ui+y^5(A'f) (8> 
Ui-y45(A^f) S U(2A^f) generated by ^®(l±r45), b = \,...,Nf^ 
(c.f. Eq. ([Toll), but leaves the discrete space-time symmetries 
C, V, and T intact, in agreement with the analogous discus- 
sion in the context of QED3 iUtIi . A nonvanishing mass m in 
contrast does not break the U(2A^f) symmetry, since 745 anti- 
commutes with 74 and 75. However, as can be read off" from 
Eqs. ( fT2l - [T4b . such a mass term is odd under parity inversion 
since 745 anticommutes with P^. Because of {745, Q) - 
and (745)^ - -745 it is even under charge conjugation; since 
[745,7";,] = and (i745)* = 1745 it is also even under time 
reversal. 

A complete basis of the 4x4 Dirac algebra is given by the 
16 matrices 



= {l' r^' "^/^v/ V2, i7,,74, 17^75, 



74,75,745 



(16) 



where we have introduced the generators of the Lorentz 
transformation of the four-component Dirac spinors cr^y :- 
j[yf2,yv]- In Eq- (fT6l l. we only count those matrices cr^v with 
fi < V. A bilinear 4^7^^ is invariant under U(2A^f) transfor- 
mations if and only if 7^ anticommutes with the generators of 
Uy^(l) and 11^5(1) while it commutes with the generators of 
Ui(l) and Uy^^(l). Obviously, this is only the case for ijryfiij/ 
and 1^7451/'. Imposing an invariance under U(2A^f) as well as 
C, and T, there is no bilinear to zeroth derivative order In 
particular, no mass term is permitted. To first order, only the 
standard kinetic term 



■kin 



(17) 



can appear Consequently, on the level of four-fermion in- 
teractions, the Thirring interaction is not the only fermionic 
4-point function in the point-like limit (i.e., with momentum 
independent couplings) which is invariant under the present 
U(2A^f) flavor symmetry and the discrete space-time symme- 
tries. In fact, the possible interactions are 



and 



(^745-^") 



(18) 



^ Note that due to our chiral conventions a nonzero expectation value (i//ifi) 
or {1^745 1/') is purely imaginary L4(]il . 



as well as the two interaction terms with non-singlet flavor 
structure 

y2 := {rif^'f - [ry^'p'f - {rrstf^'f + {ry,5^')^ , (19) 



(«A"7;.^')'+(^"^<A'') -(ri7,74/)'-(«A"i7,75«A'')' , 



(20) 



where we define (i/r°i/''')^ = ^"ij/''^^^", etc. However, the 
terms in Eqs. ( fT9] ), ( |20l ) are not independent of those in 
Eq. ( fTSl l but can be mapped onto each other by means of 
Fierz transformations. We find that the scalar/pseudoscalar 
interaction in Eq. ( fT9l ) is equal to - -Sj^ - and the 
axial/vector-type interaction in ( |20l i obeys - Sj, - 35^; see 
Appendix lAl This generalizes the discussion of |(7|] to larger 
flavor number A^f . 

To summarize: in addition to the Thirring coupling 5^, a 
second point-like linearly independent four-fermi coupling 5^ 
satisfies the symmetries of the Thirring model. In an RG anal- 
ysis, it has to be included on the same fundamental level as 
the Thirring interaction. 



III. FERMIONIC RG FLOW 

All physical information of a quantum field theory is stored 
in correlation functions which in turn can be extracted from a 
generating functional. By a Legendre transform of the latter 
one obtains the effective action F, which governs the dynamics 
of the field expectation value, taking the effects of all quan- 
tum fluctuations into account. In other words, a given the- 
ory is solved, once F is computed. Instead of integrating out 
all fluctuations at once, we can implement Wilson's idea of 
integrating out modes momentum shell by momentum shell, 
leading us to the scale dependent effective average action Fit, 
with a momentum-shell parameter k. corresponds to the 
bare action S - J d^x X. for k approaching the UV cut-off A, 
while the full quantum action F is approached for — > 0. The 
scale dependence of F^ (as a functional of only fermionic de- 
grees of freedom in our case) is governed by the Wetterich 
equation jssll 



5,F,[iA,i^] =-^Tr 



kj 



d, = k 



d_ 

dk' 



(21) 



where the trace is meant to be taken over all internal degrees 
of freedom (flavor, spinor, momentum). Here, F*'^^[0', t//] is the 
second functional derivative with respect to if/ and tp, and Ri^ 
denotes a momentum-dependent regulator function, ensuring 
that IR modes below the momentum scale k are suppressed. 
The minus sign on the right hand side of Eq. (1211 1 is due to the 
Grassmann nature of iff and if/. For reviews, see e.g. [41, 42]. 
With the Wetterich equation being an exact equation, consis- 
tent approximation schemes can be devised that allow for a 
systematic nonperturbative investigation of the given model. 
In this work, we use a simple derivative expansion of the ef- 
fective action in terms of purely fermionic degrees of freedom 
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with point-like interactions, 

(22) 

In addition to the interaction terms discussed above, we have 
included a wave function renormalization Z^. All parameters 
in the effective average action are understood to be scale de- 
pendent which is indicated by the momentum-scale index k. 
This truncation corresponds to a next-to-leading order deriva- 
tive expansion, which can consistently be extended to higher 
orders and thus defines a systematic nonperturbative approx- 
imation scheme. As discussed in the previous section, the 
truncation (l22l i represents a full basis of fermionic 4-point 
functions in the point-like limit, which are compatible with 
the present chiral and discrete symmetries. Such a point-like 
truncation can be a reasonable approximation in the chirally 
symmetric regime, as has been quantitatively confirmed for 
the zero-temperature chiral phase transition in many flavor 
QCDiH. 

Inserting Eq. (l22T i into Eq. ( 1211 1. we obtain the flow equa- 
tions (i.e., /3 functions) for the 4-fermi couplings gi^ and gi^ and 
the wave function renormalization via suitable projections 
onto the associated operators. For the explicit computations, 
we refer the reader to Appendix iB] In terms of renormalized 
fields 



we obtain ^ - 1 . For a linear cut-ofi^ which satisfies a regu- 
lator optimization criterion [43], 



and dimensionless renormalized couplings 

g = zjk-'gt, g = z^k-'-g,, 

we obtain the beta functions as 



9,8^8 — 



dt8^8- 



At 



<F) r 



INj - 1 ,2 3 

2Nf + 1 2 



2Ni 
1 



Nf 



6Nf 



(23) 
(24) 

(25) 
(26) 



Within the present truncation of point-like interactions, the 
anomalous dimension remains rji^ := -d. In Zk = 0. These 
flows involve the threshold function f^p which encodes the 
details of the regularization scheme as specified by the dimen- 
sionless regulator shape function r{q^/k^) defined by Rkiq) = 
Zk(iriq^/k\ 



d_ 

'dkJo 



d\q\ 



1 



[l+riq^k^)] 



2^2 ■ 



(27) 



For a given regulator function, this integral in the present trun- 
cation boils down to a simple number. For instance, for the 
sharp cut-off" 



I oo for q < k 
lO for q^>k^ 



(28) 




(29) 



we obtain /j^' = 2/3. By another rescaling g i-> gjr^/Atf ^ and 
g i-» g7T^/4f^p, this multiplicative regulator dependence drops 
out. 



dig = g + 



1 

2Nf^ 



2N{+1 2 
6A^f ^ ■ 



(30) 
(31) 



For A^f - 1, our result coincides with the result found via 
a perturbative RG approach |7] in the context of interacting 
electrons on the honeycomb lattice. 



IV. FIXED POINTS AND CRITICAL EXPONENTS 

From the coupling flows, it is straightforward to analyze the 
fixed-point structure in order to study possible asymptotically 
safe UV trajectories of the RG flow. A fixed point g* is defined 
by 



V/: Pi(g\,gl,...)^0, 



(32) 



with /3i = dtgi- Whereas the fixed-point values are regulator- 
scheme dependent, see above, the linearized flow in the vicin- 
ity of the fixed point is universal as quantified by the criti- 
cal exponents. More specifically, the Jacobian B,-' of the lin- 
earized flow near a fixed point. 



d,gi^BiKgj-g))- 



d/3i 



(33) 



defines the stability matrix. The associated eigenvectors v 
govern the evolution of small deviations from the fixed point 
according to d,v = Bv - -&v. The corresponding eigenvalues 
(including a minus sign) © are universal and can be associated 
with thermodynamic critical exponents if the fixed point cor- 
responds to a critical point of a 2nd order phase transition. For 
brevity, all 0's are referred to as critical exponents. The solu- 
tion to the linearized flow v oc A: ® implies that positive > 
correspond to RG relevant, i.e., infrared repulsive, directions 
and negative < correspond to RG irrelevant, i.e., infrared 
attractive, directions. 

For the present set of flow equations, a general property 
of the fixed points and their critical exponents can be proven 
iHl : the beta functions all are of the form 



A- = 8i + 8kAi8i 



BJ = 5J + 2g*Af^ (34) 



with matrices Af which are symmetric in the upper indices. 
We now see that for every interacting fixed point g* the 
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Nf = 100 Nf 
,.J. 




FIG. 1 . Positions of RG fixed points and (ir-)relevant directions for 
flavor numbers A'f = 1,2, 4, 10, 100. Arrows denote the RG flow 
towards the IR. In agreement with Eq. J35t . one relevant direction of 
every interacting fixed point points towards the GauBian fixed point 
O. 



fixed-point vector itself is an eigenvector of B, v = g* , with 
the critical exponent 0=1, 



BJg)^g\^2g*^^>g)^-g 



(35) 



Here, we have made use of the fixed point equation /3j(g*) - 0. 
We conclude that every fixed point besides the GauBian fixed 
point ^* = has at least one relevant and thus infrared re- 
pulsive direction. Each non-GauBian fixed point is therefore a 
candidate for a possible UV completion, potentially defining 
an own universality class ll44ll . For the GauBian fixed point 
g* = 0, the stability matrix is just the identity Z?,-' = 5,-', such 
that the GauBian fixed point is infrared attractive in every di- 
rection with = -1, giving rise to only trivial theories at long 
ranges. 

For fixed gj^i, the beta function yS, corresponds graphically 
to a parabola, such that we expect for our truncation exactly 
2^ = 4 (possibly complex or degenerate) solutions of the 
fixed-point equations. For any A^f e N, we find them to be 
real and non-degenerate; the explicit solutions ig*,g*) for the 
GauBian fixed point O and the three non-GauBian fixed points 
^,S,Care 



O : (0,0) 

2M- 



y2Nf- r 



(36) 
(37) 



A^f (l4 - INf + 2Nf^ + (1+ 2Nf) ^J 16 + 28Nf + Nf^) 



-5 + Mt + 2Nt^ + 4Nf'^ 



4 + 4Nf + 4Nf^- - ^Jl6 + 28Nf + Nf^, 



(38) 




Nf 



FIG. 2. Critical exponents as a function of flavor number A'f for 
the interacting fixed points S, C. The value of the corresponding 
other critical exponent is = 1 for all A'f. Note that the scale of the 
vertical axis changes at = - 1 . 



C: 



'M- (l4 - 7A^f + 2Nf^ -{1+ 2Nf) yjl6 + 2mf + N{^) 



-5 + 8M- + 2Nf^ + 4N/ 
l2Nf^ ^ 



4 + 4Nf + 4N{^ + Vl6 H- 28M + A^f^ 



(39) 



It is straightforward to derive the eigenvectors and eigenval- 
ues of the stability matrix (and thus the critical exponents and 
the RG (ir-)relevant directions) analytically; but as the general 
formulas may not provide much physical insight, we present 
the results graphically: in Fig.[Tl we plot the positions of the 
fixed points in theory space (spanned by the two couplings 
g and g) together with the corresponding eigenvectors of the 
stability matrix B,-' for various flavor numbers A^f. In Fig. |2] 
the (non-trivial) critical exponents for the interacting fixed 
points are given as a function of A^f . For any A^f , fixed point S 
has two RG relevant directions, whereas the interacting fixed 
points ^ and C have one relevant and one irrelevant direc- 
tion. As the number of relevant directions corresponds to the 
number of physical parameters to be fixed, theories emanating 
from J?[ and C are fully determined, once the initial condition 
for this relevant direction is fixed. In the sense of dimensional 
transmutation, fixing this one parameter can be viewed as fix- 
ing a total scale for the system. Therefore, theories belonging 
to these universality classes defined by ^ and C are fully pre- 
dictive, once a global scale is fixed. Theories emanating from 
S are fixed by a mass scale and one further parameter, e.g., a 
dimensionless coupling ratio, whereas O does not support an 
interacting system at long ranges. 

At this point, let us already stress that no interacting fixed 
point is on the pure Thirring axis (g = 0) for any finite A^f . In 
fact, as discussed below, we associate the Thirring universality 
class with fixed point C which approaches the pure Thirring 
coupling only in the asymptotic limit A^f — > oo. For any finite 
A^f, the renormalized UV trajectory of the Thirring model will 
have to pass through the full two-dimensional coupling plane, 
even though the long-range physics does depend only on one 
physical parameter (e.g., the value of the Thirring coupling at 
a certain scale). 
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FIG. 3. Classification of Thirring-like 4-fermi theories determined 
by tlie fixed-point positions and the corresponding RG trajectories 
(arrows denote the flow towards the IR). Red (solid) lines depict sep- 
aratrices that interpolate between fixed points and separate different 
regions. The vertical axis g = corresponds to models with a pure 
Thirring coupling. The dashed line {g = g/2) marks a generalized 
NJL model, see text. 



V. PHASE TRANSITIONS AND LONG-RANGE PHYSICS 

A technical means for the discussion of long-range phases 
are the separatrices, i.e., those RG trajectories that interpolate 
between two fixed points. They subdivide the theory space 
into separate flow regions, providing a classification which 
can potentially be related to spontaneous symmetry break- 
ing in the long-range limit. At this point, we stress that the 
fermionic truncation is not sufficient for a complete discus- 
sion of long-range physics which is expected to be dominated 
by composite bosonic degrees of freedom such as condensates 
and excitations on top of condensates. 

Let us start with a closer look at the one-flavor case A^t = 1 ■ 
The separatrices subdivide the theory space into distinct sec- 
tions defined by their IR and UV behavior; in Fig.[3l we plot 
the RG flow using Eqs. (I30b-(|3T]). A classification of RG tra- 
jectories is listed in Tab. U The IR behavior of the theories in 
the regions I, Ilia and Illb is governed by the GauBian fixed 
point O. In these regions, both couplings g and g are irrele- 
vant, leading to non-interacting theories in the IR. The regions 
Ila and IVa are characterized by an irrelevant Thirring cou- 
pling g, lim^^^o § = 0. We thus expect that the bosonic channel 
S ~ ijr"'/4sij/" becomes critical at a sufficiently large g, dynam- 
ically generating a parity breaking mass m cc {i^yn^ij/). If so, 
the fixed point governs the spontaneous breaking of parity, 
potentially being associated with a 2nd order phase transition 
defining a new universality class. 



TABLE I. Classification of all RG trajectories, c.f. Fig.[3] 



Type 


UV behavior 


IR behavior 


universality class 




limt.^o(g, g) 


limt^A(g,^) 




I 


•B 





non interacting 


I-IIa separ. 


S 






I-IIb separ. 


S 


c 




I-IIIa separ. 







non interacting 


I-IIIb separ. 


C 





non interacting 


Ila 


S 


(oo,0) 


parity breaking 


Ila-IIb separ. 


S 


S-oo 




Ila-IVa separ. 




(oo,0) 


parity breaking 


lib 


s 


C- oo 


chirality breaking 


Ilb-IVb separ. 


c 


C- oo 
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By contrast, both g and g diverge in the regions lib and 
IVb in the infrared limit. To interpret this behavior, let us (for 
arbitrary A^f > 1) rewrite the interaction terms by means of the 
Fierz theorem as 

- ~g f(r - {^"yA^f - irjsrf] 

+|[(^V45'A")'-(«A>45<A'')']}, (40) 

cf. Eqs. (fT8])-(l20]i and App.|A] In the one-flavor case A^f = L 
the last term (proportional to +g) vanishes. The interaction 
then simply is a linear combination of the Thirring interac- 
tion with coupling 2? - g and a generalized Nambu-Jona- 
Lasinio (NJL)-type 04511 interaction with coupling -g. Along 
the straight line through O and C which may govern the IR be- 
havior of the trajectories lib and IVb the NJL-type interaction 
in fact dominates. If |/|2g — g\ ~ 4.24 » L To illustrate this, 
we have also plotted in Fig.[3]the line of vanishing Thirring in- 
teraction 2g-g-Q 'm the Fierz-transformed form as a dashed 
line. This resulting NJL-type line is fairly close to the separa- 
trix through O and C (red line). We thus expect that at a large 
coupling on this separatrix the NJL-type channel eventually 
becomes critical, inducing a symmetry-broken state with 

m cc {if/i//} + 0. 

This is equivalent to saying that {xj/y^tli) + or {^'ji^') is 
expected to be favored in this state. For A^t = 1, we there- 
fore identify the fixed point C with the critical point gov- 
erning the phase transition into the chiral symmetry broken 
phase for all theories of the regions lib and IVb, in agreement 
with liTJ. This is precisely the behavior which is expected in 




FIG. 4. Fixed-point positions and separatrices in the (g, g) coupling 
plane as a function of flavor number for 1 < Nf <7. The horizontal 
slice at A'f = 1 is equivalent to Fig. [3] 



the Thirring model, hence we associate all trajectories ema- 
nating from C with UV complete fully renormalized versions 
of the 3d Thirring model. 

Provided the last term in Eq. ( |40l i proportional to +g does 
not dominate even for higher number of flavors, we may ex- 
tend the preceding discussion to larger A^f > 1. In Fig.|4] we 
show how the positions of the fixed points and the separatrices 
behave for increasing flavor number For A^f — > oo the region 
I turns into a rectangle with the vertices (§, g) - (0, 0) and 
(1, -3). The Thirring fixed point C hence moves to (0, -3) 
and becomes a pure Thirring coupling, c.f. also Fig.[T] Thus, 
the IR attractive line OC approaches the Thirring axis g - 
for larger flavor number. 

Along this line of coupling values, we expect the flow to 
be dominated by the vector channel tf/'j^jil/" which agrees pre- 
cisely with the dominant bosonic degree of freedom in a large- 
A^f analysis ll24l - l27ll . As there is no chiral symmetry breaking 
at large A^f, it is natural to expect a quantum phase transi- 
tion to occur for increasing A^f while the line OC undergoes 
a transition from the NJL-regime g ~ g/2 to the large-A^f 
regime where g = 0. Unfortunately, a more quantitative pic- 
ture of this quantum phase transition is difficult to obtain in 
the purely fermionic language. A quantitative RG analysis re- 
quires the inclusion of dynamical chiral (i.e., NJL-type) and 
vector bosonic degrees of freedom in order to study the inter- 
play of these competing orders as a function of A^f . This is left 
for future work. 

For the remainder of this section, we shall be satisfied with 
a simple estimate of the transition region. As a rough criterion, 
let us determine the flavor number where the separatrix OC is 
half way in-between the generalized NJL model on the one 
hand and the pure Thirring coupling on the other hand. This 
value of A^f follows from 



III 



\2g-~g\ 



Nf-7 + ^Jl6 + 28Nf + Nf^ 



1, 



(41) 



FIG. 5. Comparison of non-universal critical Thirring couplings 
from different methods. Solid lines display the critical couplings 
from this work taken as the g coordinate of the intersection point 
of the separatrix SC with the Thirring axis g = for the sharp cutoff 
(upper solid/brown line) or the linear regulator (lower solid/red line). 
Dotted: Monte Carlo results (magenta with circles: |30], purple with 
squares: fsT]). Dashed: DSE approaches (from right to left: green 
L27.1. cyan [28J and blue 1.29.1 ). 



critical behavior can easily provide further coiTection factors 
of order C)(l). 

It is instructive to compare our results with those obtained 
by other methods such as Monte Carlo simulations or trun- 
cated DSEs. A variety of studies have computed the mod- 
ulus of the critical Thirring coupling which is a necessarily 
required for chiral symmetry breaking (but not suflicient be- 
yond the critical flavor number). Two cautionary remarks are 
in order: first of all, these critical couplings similarly as the 
fixed-point values are not universal, such that the choice of 
the regularization scheme can have a strong quantitative in- 
fluence on the estimated values. Second, most other studies 
have defined the microscopic Thirring model exactly on the 
ThiiTing axis § = 0; in principle, the full coupling plane has to 
be considered such that the critical coupling on the axis may 
be different (larger in modulus) from a coiTesponding estimate 
directly at the fixed point. 

To circumvent the second caveat, we also consider the cou- 
pling on the Thirring axis g - with an initial value of the 
coupling g above the separatrix which interpolates between 
the fixed points S and C. Then the theory is in the region 
I or Ilia, depending on the sign of g, such that the RG flow 
drives the couplings to the free theory at GauBian fixed point 
O. Therefore, the absolute value of the g coordinate of the 
intersection point of this separatrix with the Thirring axis pro- 
vides a lower bound for the absolute value of the critical cou- 
pling gci- at which the ;^SB phase transition occurs. We com- 
pare this lower bound for different regulator functions 7?^ with 
the values of the critical coupling obtained by Monte Carlo 
simulations I30ll3ll l and different sequences of truncated DSE 
In Fig. |5] we plot this inverse coupling 
< A^f < 6. For the 



31 



and is given by A^f » 7/4. We stress that this number should 
not be viewed as a direct estimate of A^"", as the onset of true 



pc, 

approaches 11271 - 

i/ga- for varying number of flavors 
sharp cut-off as well as for the linear cut-off our results lie 
well above the values obtained by DSE approaches but below 
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the values from lattice-regularized Monte Carlo simulations. 
Note that similar to the lattice and one of the DSE results, we 
do not observe a sharp decay of the inverse critical coupling. 
This behavior would indicate a sharp growth of above a 
certain number of flavors, corresponding to a critical flavor 
number in the infinite coupling limit. Instead, we observe a 
rather smooth dependence on A^f as on the lattice for A^f < 6, 
which is compatible with our expectation that the quantum 
phase transition towards the chirally symmetric phase occurs 
because of competing large-A^f degrees of freedom and not be- 
cause of a change in the UV critical structure. 



VI. CONCLUSIONS 

In this work, we have investigated 3d relativistic fermionic 
models in a theory space, defined by chiral and a set of dis- 
crete symmetries and point-like interactions. Even though the 
construction of these models has been inspired by the uniquely 
fixed 2d Thirring model, the corresponding 3d symmetries in- 
volving a reducible 4-component Dirac spinor representation 
enlarge the minimal coupling space and give room for a larger 
fixed-point structure and thus for different microscopic real- 
izations of such fermion models. 

We have classified all point-like interactions satisfying the 
symmetry constraints and determined their RG flow in a sys- 
tematic next-to-leading order derivative expansion. The fact 
that leading-order and next-to-leading order results are iden- 
tical as the anomalous dimension remains zero can be inter- 
preted as a signature for the convergence of the expansion, as 
long as no further composite channels develop a strong RG 
flow. 

The resulting flow equations for the two independent fer- 
mionic couplings generate a fixed-point structure of one triv- 
ial GauBian and three interacting fixed points yi, S, C which 
can be classified by their critical exponents. We associate all 
RG trajectories emanating from the fixed point C with fully 
renormalized UV complete versions of the 3d Thirring model; 
as this fixed point has one RG relevant direction, fixing one 
physical scale suffices to obtain a fully IR predictive quantum 
field theory in the Thirring universality class. 

Fixed point which, incidentally, lives in an RG invariant 
sub-manifold of theory space (defined hy g - 0) is also char- 
acterized by only one RG relevant direction. From the nature 
of the scalar channel associated with this coupling direction, 
we conjecture that this fixed point can be related to a critical 
point of a 2nd order phase transition to a phase with broken 
parity. The third interacting fixed point S has two relevant 
directions; depending on the initial conditions, theories ema- 
nating from this fixed point can flow to the Thirring phase as 
well as to the parity-broken or symmetric phases. 

Unfortunately, our purely fermionic RG analysis does so 
far not permit us to reliably run towards or into the symmetry- 
broken phases. Such a quantitative description is required in 
order to analyze the true IR behavior of the Thirring phase. 
From the structure of the fermionic flow, in particular, from 
the behavior of the separatrix OC, we conjecture that the long- 
range dynamics is characterized by a competition between 



NJL-type chiral condensation channels on the one hand and 
large-A^f-type vector bosons on the other hand. As the vector- 
boson fluctuations generically inhibit chiral symmetry break- 
ing, we expect the occurrence of a quantum phase transition of 
the Thirring model at a critical flavor number N", separating 
a broken phase for small A^f from a disordered phase for large 
A^f. Our very rough estimate of the transition region in any 
case is compatible with the findings from lattice simulations 
indicating that A^° ^ 6.6 Isill. 

Our results on the position of the Thirring fixed point C 
being actually away from the pure Thirring axis g = for 
any finite A^f provokes an important comment: both lattice 
simulations as well as DSE studies build on a microscopic 
definition of the 3d Thirring model which is fixed only with 
the Thirring coupling, i.e., by a pure bare Thirring-like ac- 
tion. Our fixed-point results indeed provide a fundamental 
justification for this, as the Thirring fixed point indeed is char- 
acterized by only one relevant direction. As long as the mi- 
croscopic actions chosen in other formulations are in a suffi- 
ciently attractive domain of the Thirring fixed point, univer- 
sality guarantees that the long-range physics is indeed purely 
governed by the Thirring fixed point. We expect this to hold 
also for the determination of the critical flavor number (there 
is no universality for quantities such as the scheme-dependent 
critical coupling, see Fig.|5]l. Nevertheless, one caveat should 
be emphasized: this conclusion about universality only holds, 
as long as the microscopic bare actions indeed are in the at- 
tractive domain of the Thirring fixed point. For instance, if by 
accident a lattice formulation turned out to be influenced by 
fixed point S, the simulation would simply describe a model 
different from the Thirring model potentially exhibiting a dif- 
ferent quantum phase transition as a function of A^". As the 
fixed-point positions are not universal, our results are unfor- 
tunately not directly transferable to the lattice theory space. 
However, provided that the fixed-point structure is qualita- 
tively similar, our results can be taken as a support for the 
implicit assumption that the lattice simulations have indeed 
been performed in the real Thirring universality class. 

Given the importance of the quantitative value of the crit- 
ical flavor number A^" of the Thirring model in the light of 
condensed-matter applications, a natural next step of our stud- 
ies will be the inclusion of composite degrees of freedom 
in order to study the competition among the various bosonic 
channels. Within the functional RG, this can conveniently be 
investigated by means of partial or dynamical bosonization 
[46] . As the problem of competing order is a paradigmatic one 
in (quantum) critical phenomena and statistical physics |47], 
we consider the relativistic 3d Thirring model as an ideal test 
case. 
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Appendix A: Fierz transformations 

Let (7^,7^) '■- Tr(y^7^) define an inner product on the 
space of 4x4 Dirac matrices. The 16 matrices {y^} in Eq. ( fT6] l 
are orthogonal with respect to this product, 



Tr(/r*) = 46' 



AB 



(Al) 



and thus define a complete basis of the 4x4 Dirac matrices 
(which generically have 16 independent matrix elements). 



16 J 



(A2) 



The Fierz transformations are straightforwardly obtained by 
multiplying the completeness relation (IA2l l by each of the 4- 
fermi terms \j/1J,yA^^')k^'j{yA^'')i and, where appropriate, de- 
composing multiple products of Dirac matrices into basis ele- 
ments jA ■ We find 

16 
B=l 
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(Cab) = 



-1 


-1 


-1 


-1 


-1 


-1 


-1 


-r 


-3 


1 


3 


1 


-1 


-1 


-3 


3 


-1 


1 


-1 


-1 


1 


-1 


1 


1 


-3 


1 


-3 


1 


1 


1 


-3 


-3 


-3 


-1 


3 


1 


1 


-1 


3 


-3 


-3 


-1 


-3 


1 


-1 


1 


3 


3 


-1 


-1 


1 


-1 


1 


1 


-1 


1 


-1 


1 


1 


-1 


-1 


1 


1 


-h 



(A4) 



tween the invariant 4-fermi interactions in Eqs. ( fTsT i, ( fT9] ) and 
(Unil. 



Appendix B: Derivation of flow equations 

Consider the right-hand side of the flow equation ( l2Tl i. By 
introducing the derivative d,, acting per definition only on the 
f dependence of the regulator Rf., we may expand 



5,in(r[^ + /?, + Arf) 
5, In (rg + + a, 

' AT? " 



(Bl) 



Ar 



<2) \ 



<2) 



(B2) 



with AF^ containing the field-dependent parts of F^ , and 
AF^ Q containing the field-independent (propagator) part, such 



that F' 



<2) 



FJ,Q + AFj. . The Hessian of the average effective 



action in terms of Fourier transformed fields ij/iq) = ijfq and 
4i{q) = ipc, is given by 



6 S 6 6 

6 6 6 y 6 



(B3) 



and (ja) = (1, 74, ir/ir4, iT/jTs, 745, Ts)^- With these 
preliminaries one can simply read off the Fierz relations be- 



The corresponding fluctuation matrix results in 



AFf(p,g) 



gk_ 
Nf 



+ term with (gt, 745) <^ (gk, Jp) 



L {(<^.^45)^(r45«A*-,.p/ - {^p.yA5^p-,.p)jlS"''] 

-/„(r45r„)(r45^;;-,-„f 



(B4) 



with s J and for the propagator 



ZuqHl+riq^lk'-)) 



q 

f 0, 



withRkiq)^Zkqr(qyk^). 

Anomalous dimension. The flow of the wave function 
renormalization is obtained by a suitable projection of the 



r 



Wetterich equation (1211 1. 



(B5) 5^z, = -Tr^r, 



- f — 

dp'p Jq' Sifrl 



Tr 



d,Rk 



\ rf * + Rk) 



sr. 



p ^q- -Tp, 

Using the expansion ( 

(2) 

in AF^ survives the 
the flow of the wave function renormalization in this order of 



i/(=i//=0 ) 

(B6) 

, we observe that only the term linear 

(2) 

in AF^ survives the projection. Since this term is traceless. 



10 



the expansion vanishes, 

T]k = -d,lnZk = Q. (B7) 

This line of argument is reminiscent to a similar observa- 
tion of the vanishing of the anomalous dimension in scalar 
0{N) models in the symmetric regime at next-to-leading or- 



der derivative expansion. 

4-fermi couplings. Obviously only the term in Eq. ( IB2I) 
being quadratic in AF^'^' leads to 4-fermi terms on the 
right-hand side of the Wetterich equation ( 1211 1. contributing 
to the flow of the 4-fermi couplings. Evaluating this term for 
constant fields and taking the trace over flavor, spinor, and 
momentum degrees of freedom, we infer 



^ k,0 + 



2Q 

Nfn- 



d\q\~d, 



Zl [1 + r{q^lk^)f 



2Ni - 1 =2 



2Ni 



-8k 



2Ni 



gkgk 



Nf 



^2 I (0.^45^)2 + 



2Nf+ 1_21,- ,2 

gkgk — g^^^i I ('Ar/i'A) 



(B8) 



with Q being the spacetime volume. The desired beta 
functions for gi; and gt^, respectively, are simply given by 
2Nf/Q. times the coefficient of the corresponding 4-fermi term 
in Eq. ( IB8b . For s 1 and by taking the regulator-dependent 



constant i\ from Eq. ( |27] | into account, we end up with the 
flows of the dimensionless couplings g and g (Eq. jTM ) as 
displayed in Eqs. (|25l l and i 
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